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The microscopic spectrum of the QCD Dirac operator is shown to obey random 
matrix model statistics in the bulk region of the spectrum close to the origin using 
finite- volume partition functions. 



1 Introduction 

The spectral density p(A) of the Dirac operator in QCD contains interesting 
informations as it is for example directly proportional to the chiral condensate 
£ at the origin through the Banks-Casher relation u. While the full spectrum 
is only accessible numerically in QCD on the lattice many analytic results for 
parts of the spectrum have been obtained during the past years using random 
matrix theory (RMT) , chiral perturbation theory and finite- volume partition 
functions (for a recent review see □) . Within these results two different regimes 
have been investigated: (i) unsealed macroscopic correlations and (ii) micro- 
scopic correlations between eigenvalues rescaled by the mean spectral density. 
In the region (ii) one furthermore has to distinguish between scaling at the 
origin and in the bulk of the spectrum. In the macroscopic regime (i) the slope 
of the spectral density at the origin 

P'(A = 0) = (!) 
Hy ' 1671-2^4 pN f ' y ' 

has been calculated by Smilga and Stern i for fundamental SU(N C > 3) 
fermions {(3 = 2) and very recently by Toublan and Verbaarschot □ for funda- 
mental SU(2) fermions (/3 = 1) and adjoint fermions ((3 — 4) with any gauge 
group. Here F„ denotes the pion decay constant. 

Coming to the rescaled correlations (ii) at the origin all higher order cor- 
relation functions of scaled eigenvalues x = A/EV are known. They are 
given analytically as functions of the chiral condensate E and scaled sea-quark 



masses fj,f = m//EV (for_ details see □). This description holds in the limit 
of Leutwyler and Smilga U 1/Aqcd <C V 1 ^ -C l/m^, where is the pion 
mass, up to a certain scale Xth ^V 1 / 2 called Thouless energy til in analogy to 
disordered systems. Within this region the analytic description provides an 
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Figure 1. The unsealed spectral density plotted schematically (compare Ref. 8 ). 



exact nonperturbative limit of QCD, dealing, however, with an unphysically 
heavy pion that does not fit into the finite- volume of the system. 
Away from the origin much less is known about correlations on the micro- 
scopic scale. Up to now it has been found only empirically that also there 
correlations follow a (non-chiral) RMT description LrEI The aim of the results 
presented here is to understand analytically why this is possible. Therefore, 
the ccojrelation functions are expressed in terms of finite- volume partition func- 
tions □ and explored in the bulk region but still in the vicinity of the origin. 
The region of validity is therefore given by the one of finite-volume partition 
functions. It does not hold all the way up to the unphysical turn-over in the 
spectrum (dotted line in Fig. [I]) which is due to cut-off and finite-size effects 
on the lattice. This intermediate region remains unexplained so far. 

2 Correlation functions from finite-volume partition functions 

Before taking the bulk limit we have to recall how the correlators can be ex- 
pressed using the Leutwyler and Smilga partition functions in the microscopic 
limit at the origin. The generation function for the correlators is the kernel of 
orthogonal polynomials. For gauge group n 5'[/( N c > 3) with fermions in the 
fundamental representation it is given by @ 

K s (x,y;{^}) = c^y^ + ^)(z/ + g^^fr^ iy) (2) 
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2v ({/x/}) = — ^yTj — . -An = Mi ^+j-i(/A t ) , (3) 

where zi Nf ^ is the partition function of Nt flavors with rescaled masses fif, 
v zeromodes and A being the Vandermonde determinant. The rescaled cor- 
relation functions can then be obtained as 

p { f f ' v \x u ...,x k ;{fj, f }) = det Ks{x a ,x b ;{n f }) . (4) 

l<a,b<k 

This description is equivalent to a unitary chiral RMT description (/3 = 2). 
Since we want to go to the bulk of the spectrum we should expect that the 
chiral properties no longer play a role for the correlations. We will find that 
this is indeed the case even in the vicinity of the origin. It is therefore in- 
structive to compare to 3-dimensional QC-D with flavor symmetry breaking, 
which is described by non-chiral RMT £3. The corresponding kernel can be 
written in the same way as in Eq. (0) a after dropping the prefactor c^/xy, 



where the product now runs over Nf/2 flavors occurring in pairs ±/x/. The 
corresponding partition function reads 

*•*•«*» - « ($$ } , = («rv.. (5) 



3 The bulk limit 

In RMT it is known for massless flavors how to "invert" the origin scaling 
limit where V — > oo, A — > and x = X/T,V is kept fixed: 

lim K OIlEin {x,y) = K bulk (x,y) = d — i ^ . (6) 



x, y — > oo 
x-y = 0(l) 



x-y 



It takes us from the Bessel-kernel at the origin (Eq. (||) for m/ = 0) to the 
sine-kernel in the bulk. We will see in the following that the same result holds 
including the masses. In RMT the translational invariance of the spectrum 
then implies that the sine-kernel holds everywhere in the bulk. Since we 
cannot make such a strong statement about the QCD Dirac spectrum we will 
be restricted to the domain of validity of finite-volume partition functions. 
In taking the bulk limit of the full kernel Eq. (^) we also have to scale back 
all masses [if — > oo at the same time since they would otherwise drop out 
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trivially. Using the asymptotic expansion of Bessel functions we obtain 

hm K s {x,y;{Lif}) = c [[ — - .(7) 

x,y,fi f -> oo f=1 J (x 2 + n)){y 2 x y 

difforcncoa O(l) 

Expanding the prefactor around one of the arguments it is unity up to higher 
orders and thus we find the same result as in the massless case Eq. (Q) . Taking 
the same limit of the non-chiral theory Eqs. @,(|5|) we obtain precisely the 
same result and thus find, that the chiral properties are lost in the bulk even 
in the vicinity of the origin. Another consequence is that in the bulk any mass 
scale giYen by the sea-quarks drops out on the microscopic scale. 
In Ref. LJ it has been stated that the Thouless energy scale at the origin is 
larger or equal to the one in the bulk. Looking at the bulk but close to the 
origin it is clear that if we take Ai < At/i and A2 > Xth we are outside the 
domain of validity although it may still be |Ai — A2I < Xth- 
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